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We present a new formulation for the evaluation of the primordial spectrum of curvature pertur- 
bations generated during inflation, using the fact that the Wronskian of the scalar field perturbation 
equation is constant. In the literature, there are many works on the same issue focusing on a few 
specific aspects or effects. Here we deal with the general multi-component scalar field, and show 
that our new formalism gives a method to evaluate the final amplitude of the curvature perturbation 
systematically and economically. The advantage of the new method is that one only has to solve a 
single mode of the scalar field perturbation equation backward in time from the end of inflation to 
the stage at which the perturbation is within the Hubble horizon, at which the initial values of the 
scalar field perturbations are given. We also clarify the relation of the new method with the new 
8N formalism recently developed in Ref. 

PACS numbers: 



I. INTRODUCTION 



Inflation provides an elegant mechanism to solve shortcomings of the standard Big Bang Model, for 
example, the flatness and horizon problems. Moreover, in the inflationary universe, curvature perturbations, 
which seed the structure formation of the universe, are generated from vacuum fluctuations of the scalar 
field. One can test various models of inflation, by comparing the theoretical prediction for the spectrum of 
curvature perturbations with observations. 

In single-field slow-roll inflation models, curvature perturbations evaluated on comoving hypersurfaccs, 
1Z C , arc constant on super-horizon scales, and their spectrum is given by 

/ tt2 \ 2 

VnM~ (r-i) , (1-1) 

aH=k 

where H is the Hubble parameter, <fi is the derivative of the inflaton </> with respect to the cosmological 
time, and the subscript u aH = fc" means that the expression is to be evaluated at the horizon crossing 
time, t = tk- This equation indicates that the spectrum of the primordial curvature perturbations is almost 
scale-invariant in slow-roll inflation models composed of a single scalar field. 

However, in making realistic models of inflation based on supersymmctry and supergravity theories, it 
seems more natural to consider models composed of a multi-component scalar field, in which the slow-roll 
conditions are also possibly violated 0. Obviously, Eq. is not sufficient to evaluate the spectrum of 

the curvature perturbations generated in such models. In order to identify necessary conditions for viable 
models of inflation among various proposed possibilities, any simple and sufficiently accurate formula for the 
spectrum is useful, especially if the formula is applicable to a wide class of models. 

There are several works aiming at generalizingthe slow-roll formula in the case of a single scalar field. 
For example, Leach, Sasaki, Wands and Liddle [3j have analyzed the evolution of curvature perturbations on 
super horizon scales (the long- wavelength approximation). They have found that the curvature perturbations 
can be enhanced when z = acj)/ H decreases after the horizon crossing time it. This situation can happen if 
the slow-roll conditions are violated for t > tk- On the other hand, Stewart |j| has generalized the standard 
formula removing the extra assumption that slow-roll parameters are approximately scale invariant, 



and this generalization is called the general slow-roll approximation. Under the conditions which lead to 
the standard slow-roll formula the running of spectrum is automatically suppressed as \dnn c /d\nk\ <C 

\nji c — 1|, while under the general slow-roll conditions, \dnj3 c /dhik\ can be as large as \n-ji a — 1|. These 
conditions are consistent with recent observational data la] ■ 

For multi-component inflation, evaluating the curvature pert urbations is a more difficult issue in compar- 
ison with single-component inflation [1 II IE II El El El 111 El El El El E3 Generally, in the case of 
a D-componcnt field, the evolution equations for the perturbations arc D coupled second order differential 
equations. It is a heavy task to solve the full set of equations. Moreover, for multi-component inflation, 
the curvature perturbations on comoving slices, 1Z C , still varies in time on super- horizon scales until the 
neighboring trajectories of the background homogeneous universe in the phase space practically converge to 
a single one. It is after the convergence of the trajectories that TZ C becomes constant in time |l7l Il8|. 

However, in the estimation of the spectrum, what we need is only the final value of the curvature pertur- 
bation lZ c (jif). Here, 77 is a conformal time given by 77 = J dt/a, and 77 f is taken sufficiently late after the 
convergence of the background trajectories. So, if we can identify the part of perturbations that contributes 
to TZ c {r]f), we may solve only that part to obtain its value without solving the full set of perturbation equa- 
tions. The SN formalism 0, 0| was based on this idea, and we have developed a new SN formalism 0. 
This formulation was developed based on the fact that the perturbation of e— folding number stays roughly 
constant in the flat slicing on super-horizon scales. We derived a perturbation equation for a single projected 
component of scalar field corresponding to the perturbation of e— folding number. The derivation of the new 
SN formalism is rather intuitive in this sense, which helps to make its physical meaning rather transparent. 
As a drawback, however, the systematical derivation of the formula is not so simple. Moreover, in Ref. 
we only applied our new formalism to cases without the possible super-horizon contributions pointed out in 
Ref. [3j , postponing this generalized application to Ref. |l9j . Here we propose an alternative slightly different 
framework to evaluate the evolution of perturbations of a multi-component scalar field. The method we will 
discuss in this paper takes full advantage of the fact that the Wronskian stays constant. 

This paper is organized as follows. In section [HJ we briefly explain the perturbation equations for multi- 
component scalar field in the flat slicing. In section ITTT1 we introduce a Wronskian between two solutions of 
perturbation equations, and explain how we can use it in evaluating the curvature perturbation at a later 
time 77 = rjf. As far as we are concerned with the evolution during the phase dominated by the scalar 
field, what we have to do turns out to be just solving a single mode backward in time with appropriate 
boundary conditions imposed at 77 = 77 f. In section llVI we use the Wronskian to derive more explicit 
analytic expressions for the spectrum of perturbations. We consider an extension of the general slow-roll 
formula developed in Ref. 0, combining it with the long wavelength approximation d, which suites 
for describing the late time evolution |l9j . We assume that the general slow-roll approximation is valid 
until a time 77* later than the horizon crossing time 77^, and we use the long- wavelength approximation for 
the succeeding evolution until the end of inflation, 77^. To combine the two solutions obtained by using 
different approximation schemes in the matching region where both approximations are valid, we find that 
the constancy of the Wronskian is again very useful. The derived new analytic formula applies for a wider 
class of multi-component inflation models. In section we reduce the obtained formula to the single-field 
case. In this case, the formula for the spectrum can be expressed more explicitly. Using our new formulation, 
it is also easy to see that the amplitude of perturbation stays regular even if the time derivative of </> vanishes, 
as was discussed in Ref. |20j . In section I VII we present a brief summary. 



II. PERTURBATION EQUATIONS FOR MULTI-COMPONENT SCALAR FIELD 

We consider a C-componcnt scalar field whose action is given by 



d x\f—g 



\g^S aP d^ a d v ^ + V{4>) 



(2.1) 



where, for simplicity, we have chosen to work with the flat metric 5 a p for the scalar field space, because 
the equations will be more involved for the general field space metric and hence the essence of our new 
formulation may be obscured, though the extension to the general metric case is straightforward. The 



2 



background equations are 

m 2 = + a 2 V{4>) , (2.2) 

4> a " + 2H<t> a ' + a 2 S a0 V^ = , (2.3) 

where H. = a' /a, the prime represents a derivative with respect to the conformal time rj, and V^a = dV/d(f) a . 
We consider only scalar type perturbations. Then the perturbed metric is given by |2ll| 

ds 2 = a 2 {-(1 + 2AY)drf - 2BY j dT]dx j + [(1 + 2H L Y)8 lJ + 2H T Y ij ] dx l dx 3 } , (2.4) 

where Y is the spatial scalar harmonic function with the eigenvalue /c 2 , Yj = —k VjY, and Yy = 
(fc- a ViVj + |<%) Y. 
We decompose the scalar field to the background and perturbations as 

4> = <j>(n) + 8<t>(af)Y , (2.5) 

and define a gauge invariant variable 

S^f = H a - , (2.6) 

At 

which represents the scalar field perturbation on the flat slicing. Here 

U = H L + ^H T , (2.7) 

represents the intrinsic curvature perturbation of the constant time hypersurfaces. The evolution equation 
for S(j)p is known to be given by [ljj 

{d 2 - (H 2 +H >)+k 2 }^ = ^ ' S J0 - a 2 S^V^ | ^ , (2.8) 

where 

ip a = a5<f>%. (2.9) 



III. USE OF WRONSKIAN 



When we calculate the curvature perturbation spectrum, we only need the final value of the curvature 
perturbation, 7£ c (??/). We do not need to know the evolution of all components of the multi-component field. 
In other words, if we can identify the part of the perturbations that contributes to IZdrjf), we may be able 
to solve only the necessary part without obtaining a full set of solutions of perturbation equations. 

The equation (|2.8|l is an equation which formally takes the form of 

[d 2 + Q( V ) + k 2 ]^ = P a p ( v )^ 7 (3.1) 

with P a a symmetric with respect to the indices a and f3. The Wronskian 

W(n, ip) = ri ■ ip - n ■ ip', (3.2) 

introduced in a usual manner is constant if n and (f are solutions of Eq. (|3.1|l . Here we used vector notation 
n and <p to represent n a and ip a for brevity and n ■ cp = 5 a pn a LpP . 

As is known well, when the universe is dominated by a single matter component whose physical state is 
solely characterized by a single parameter, say, the energy density as in the case of a slowly rolling scalar field 
or a single component perfect fluid, the curvature perturbation on comoving hypersurfaces 1Z C is constant in 
time on super-horizon scales. If the background trajectories in the phase space of the multi-component scalar 
field well converge to a single one at a time, rfcon, before the scalar field dominant phase ends, the system 
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is effectively equivalent to a single scalar field for r\ > rj con . In such cases the problem is relatively easy to 
formulate. In comoving gauge, S a ^ a ' S^ = by definition. Contracting Eq. (|2.6[) with <j)' a = Sapcj) 13 , we 



find 



a0' 

Hence, by defining n so that it satisfies the boundary conditions 

H<t>' 



K C = (3.3) 



A.' 2 ' 

a<p 



n = 0, (3.4) 



at rj = ?7 con , the above Wronskian agrees with 1Z C : 

n c ( Vcon ) = W(ri,<p). (3.5) 

Since the Wronskian is constant, it can be evaluated at any time. Therefore, once n is solved backward 
in time until the wavelength of the mode is well inside the horizon scale with the appropriate boundary 
conditions, Eq. I|3.4[) . we do not have to evolve ip at all in the forward direction. This will be advantageous 
in the case of a multi-component field. In order to solve (p in the forward direction, we need to solve all 
2D independent modes, where D is the number of species. In contrast, to evolve n backward in time, we 
only need to solve a single mode. Hence, the advantage of solving n instead of (p becom es p rominent as 
the number of species increases. An extreme example is the models with extra-dimensions |22j. Since D is 
infinite in the brane world models, solving for all ip is impossible. Note that the boundary condition l|3.4J) is 
given solely in terms of the background quantities. This implies that n is uniquely determined for arbitrary 
scalar field perturbations. 

In general multi-component models, it may not be a good assumption to assume that the trajectories 
converge to a single one before the scalar dominant phase terminates. Once the convergence of trajectories 
occurs, the universe undergoes a universal evolution everywhere, and the only remaining modes of scalar-type 
perturbations are adiabatic ones |l7| . If we do not consider the case with isocurvature perturbations that 
persist until the present epoch, we can assume that the convergence of trajectories occurs at a certain time, 
but it need not happen during the scalar dominant phase. In multi-component models, the convergence of 
trajectories may occur in general after complete reheating 0,01- In such cases it is at this stage that 1Z C 
becomes constant in time on super-horizon scales. 

Even when the trajectories in phase space have not converged, the SN formalism [l7j simplifies the 
evaluation of the evolution of perturbations in the long wavelength limit. In the rest of this section the long 
wavelength limit also contains the meaning that we neglect the shear of the constant time hypersurfaces, 
which rapidly decays on super-horizon scales [Isj . Let us take rj i to be a certain time soon after the relevant 
scale have become greater than the horizon scale and rj f to be a certain time after the complete convergence 
of trajectories. The SN formalism is based on the fact that the final value of the curvature perturbation 
IZcivt) i s given by SAf(j)t, T]{), the perturbation in the e-folding number between the initial flat hypcrsurfacc 
at rj = rji and the final comoving hypcrsurfacc at rj = j]f in the long wavelength limit. 

The evolution of 8N~ in the long wavelength limit is expected to be described by the dynamics of the 
homogeneous universe in the following sense |l7l lljj . Wc define N by the e-folding number necessary 
to reach the point in the phase space corresponding to rjf by solving the evolution of the homogeneous 
universe. If we choose rji to be in the scalar dominant phase, N will be a function of and 0'. Then, in 
the super-horizon scales, the perturbation in N is given by 



SN 



dN dN (d<t>\ 

50 ' d<Pp+ d(dcf>/dr) \dF) 

dN /<p\ dN 
90 ' W + dfi 



a I 



(3.6) 



Here, the subscript F denotes a quantity evaluated on flat hypersurfaces and the subscript r/i denotes a 
quantity evaluated at 77 = r]i. The perturbed proper time r satisfies dr = a(l + A)dr), and the perturbation 
of the lapse function on flat hypersurfaces Ap is given by 
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Thus we obtain 



SN 



W{n B ,v) 



(3.8) 



with 



' B ~ a d<t> 



1 dN 



+ 



H dN 



+ 



1 (dN 





1 dN 



(3.9) 



2aH \d<p' 



Our expectation is that, SN is identical to SAf in the long wavelength limit in general. In fact, it was shown 
in Ref. [18j that this is indeed the case if 77 f also lies in the scalar dominant phase. We plan to discuss this 
issue in the more general cases in our future publication . Here in this paper we simply assume SAf ~ SN 
in the long wavelength limit. Notice that the meaning of this approximate equality is an identity relation 
which holds for any <p and <p' . Hence, W(n, if) « W(riB,<p) is also to be understood in the same way. 
Therefore we conclude ns ~ n. In fact, it has been shown in our previous paper that ns defined above 
is a solution of Eq. (|2.8|l in the long wavelength limit. Thus, in the case that the convergence of trajectories 
occurs after the scalar dominant phase ends, the boundary conditions for n are given by n = 713(771) and 
n' = n' B {rji) instead of Eq. I|3.4|) . 

It would be worth mentioning that n used here is a full decaying mode solution obtained without assuming 
the long wavelength limit, while in Ref. Q] we used a super-horizon decaying mode solution, which is a 
solution for Eq. (|2.8|) with k 2 — ► 0. They are identical only in the long wavelength limit. Here we make full 
advantage of the constancy of the Wronskian W{n, ip). In contrast, in Ref. we discussed the evolution of 
the super-horizon Wronskian, which is defined by W(n, ip) with n replaced by the super-horizon decaying 
mode solution. 



Even if we use the Wronskian method, we still need to solve n backward for each k. However, we can 
avoid solving n for each value of k in the long wavelength limit. In this limit, we can find the solution for n 
in the series expansion with respect to k 2 . Of course, this approximation does not hold when the wavelength 
becomes shorter than the horizon scale. For the evolution before the horizon crossing, one may use the 
slow-roll approximation or the more general formula called general slow-roll developed by Stewart 4]. By 
using these approximations, we can evaluate the evolution of mode functions to a large extent analytically. 

However, in the original slow-roll or general slow-roll formula 0,0, one assumes no contribution on the 
super-horizon scales 1 , i.e., in these formulae, we cannot take into account the super-horizon effects pointed 
out in Ref. 

Then, in order to take account of the super-horizon effects, one may naturally think of matching two 
approximation at an appropriate time 77 = 77* where the mode is already well outside the horizon scale but 
the slow-roll or the general slow-roll conditions arc still maintained. However, the matching does not seem 
to be so trivial since the solutions in the long wavelength approximation look quite different from those in 
the slow-roll or the general slow-roll approximation. 

As an application of the Wronskian method, we consider this matching problem. The goal of this section 
is to obtain an improved general slow-roll formula. Since the Wronskian W is constant in time, it is allowed 
to be evaluated at 77 = 77*. If we know both ip and n at 77 = 77*, what we have to do is simply to compute 
W(n, ip) there. Then we can avoid the messy computation for explicitly matching solutions obtained in the 
two different approximation schemes term by term. To obtain cp and n at 77 = 77* , we evolve ip in the forward 
time direction by using the general slow-roll approximation, while n in the backward direction by using the 
long wavelength approximation. 



This is roughly correct, but not exactly correct. More precisely, it assumes general slow-roll approximation is valid until 
super horizon effects disappear. 



IV. EXTENDED GENERAL SLOW-ROLL FORMULA 
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A. General slow-roll expansion 



In the single field case, under the general slow-roll condition, slow-roll parameters, 5 = 4>/H(j>, and S/H 
can be the same order but both are assumed to be small. Here we used the dot as differentiation with respect 
to the cosmological time t (dt = adrf) and H = a/ a. In this case P, which is P a p(r]) in matrix notation, 
can be shown to be small. Here we mean that max(P a g) is small by "the general slow- roll condition" as a 
simple generalization to the multi-component case. In general, not all components of the multi-component 
scalar field are necessarily nearly massless. Hence, the assumption that max(P a ^) is small is not satisfactory. 
We plan to propose an improved treatment in our future publication. We write the solution solved in the 
forward direction in an expansion with respect to P as 



<P = A V(o) + A V(i) + ••• , 
where A<p^ is the correction of 0(P P ). We also introduce variables 



(4.1) 



3=0 

They satisfy the equation 

[91 + Q(V) + k 2 -P] <p (p) = -PA<p (p) . 
Using the Green's function method, the first order correction with respect to P is obtained as 



Av(i) = n 



oo J— oo 



where n = W(uq,Uq) , and uo(rf) is a solution of 



(dl + Q{rf) + k 2 )u = 0. 



(4.2) 



(4.3) 



(4.4) 



(4.5) 



B. k 2 expansion 



The backward solution obtained by setting the final condition is obtained as an expansion with respect to 



k 2 as 



n = An {0) + An {1) + ■ ■ ■ , 
where An^ is 0(k 2q ). We also introduce 



(4.6) 



n 



(?) = 



J2An 



i=0 



They satisfy 



{d 2 + Q(n) + k 2 - P)n (q) = k 2 An {q) . 



(4.7) 



(4. 



C. Evaluation of Wronskian 



To the first order in P and to the g-th order in fc 2 , we have 

TZ[l\ = W(nM )V u)\ v .=nto'-<p w -nto-<p w ' 



W(n (9) ,¥» (0 ))k+n 



W(nM, uo)\ v . I d V 'u* (r,')P{r,')<p (0) ( V ') 



G 



W(n^,u* )\ v , / dr,'u (r)')P( V ')cp i0) (r,') 



(4.9) 



It is interesting to see how the above formula is systematically cast into the form appropriate for comparison 
with the standard slow-roll formula, The above formula still contains an ambiguous choice of the 

matching time 77*. Although the expression is independent of the choice of 77*, we cannot simply replace it 
with i]k, which is the time of the horizon crossing since the long wavelength approximation is not valid at 
?/ = r/fe. The appropriate expression is 



K 



(?) = 
(i) - 



W{n,ip (0) ) 



(«)i 



'//, 



Il{W(n,u ) / d V 'P(r,') W)¥>(o)(V)-%' -Vk)<(v')<P (0) (v') 



W(n,u* ) I tfPtf) u (v')V(o)(v')-0(v'-m)u o (v')<P(o)(v') 



(«) 



(4.10) 

Here means the truncation of the function / at 0(k 2q ) 1 while the function with tilde / means that 
we keep the original form of the function / without expanding it in powers of k 2 , namely, f^ q ' = f. This 
expression still contains 77* but the integrand takes the form of / — #(77 — rjk)f^- Therefore all the terms up 
to 0(k 2q ) cancel, and hence the integrand shows improved fall off at 77 — > 0. As a result, even if fall off of P 
in the limit 77 — > is not very fast, we may be allowed to set 77* to 0. In this sense, the independence from 
the choice of 77* is more manifest in Eq. H4.1l)fl • 

Now we will show that Ea. (|4.1U|) is equivalent to Ea. l|4.9|l . From Eqs. I|3.1fl and l|4.5|l . we can easily get 



[dr,W(n : u )} {q) =P[nu ] 



(?) 



(4.11) 



Thus, we can show that 

n • Ptp {0) 



u Vk'^(l) — 



(<y) 

Ik 



[n'uo — nuo'] P 



[nu*] P 



d^'Ptf) u (v')V>(o)(v') - Otf - Vk)u (v')<P( 0) (v') 

u 0f(0) 



[ti'Uq — UUq] P 



(4.12) 



J Vk 



The first and the third lines in the curly brackets are second order in P. The second and the fourth terms 
inside the curly brackets are combined to give 



TlW(u* ,u ) n P<p 



(0) 



u/) 



n ■ Pep 



(0) 



(?) 

Vk 



which cancels the first term on the right hand side, and hence d Vk TiM} = 0(P 2 ). 



Thus 7^ i s independent of the choice of 77^ to the first order in P. Setting rjk — > ??* in TZf, we find 



(?) 



W( V(1) ,n^%. =n[ q) +0(P 2 ,k 2 <>+ 2 ) 



(4.13) 



r 



D. Extension to the higher order in general slow-roll expansion 



The extension to the higher order in P is rather straight forward in our formulation. This clearly demon- 
strates the advantage of our systematic derivation. We introduce notations 



n[ v ] = n\w(n,u ) r dr,'P(r,')[uUv')v>W)-0(v' -VkKWMv')] 



(4.14) 



w(n, u*) / d V 'P(v') Mv'Mv') - o(v' - m)Mv'Mv')] 



(4.15) 



In the similar way as we did in the above Eqs. (|4. 1 1|> and l|4.12[) . we obtain 

■*i?fe^ /c (0) 



= [n-P V(0) ]<*>, 



(<?) 



■[n-PtpW + {nPA<pl<p} 



sn[v>] 

6(p 



(g) 



(4.16) 



where 



uS(f?) / <VP(r?') [Mv')v>W) - 0(v' - v)Mv'Mv')] 



(4.17) 

As a natural extension of Eq. I|4.10|l . we expect that 1Z including higher order corrections up to 0(P P , k 2q ) 
will be given by 



J=0 



(4.18) 



with 



Awjg=W[ V(0) ]W, 



(o)J 



(9) 



A£[A£[ V(0) ]] 
In fact, using Eqs. I|4.15[) and (|4.16() . we can show that 



(9) 



(4.19) 



,(</) 



-(o) 

-[n-PA<p[<p {0) ]]$ 



(0)J 



'(0)J 



(9) 



nPAip 



Atp[cp 



(0)1 



(0) _ 



(9) 

+ • 



0(P p+1 ), 



(4.20) 



where we used the fact that the arguments of 7Z, such as (pr Q \ and Aip[<p, \], are independent of Eq. I|4.20|l 



shows that T^jpj is independent of the choice of rjk to the p-th order in P. Setting r]k — > ?7* in , we find 



,(<?) 



(4.21) 
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E. Power spectrum 



The perturbation is to be treated as a quantum operator, and it will be expanded as 

tpfa = u a +u a 1 . 
Creation and annihilation operators have expectation value 

(a a a^) = 5 al3 . 



(4.22) 
(4.23) 



The expectation values of the other combinations of two creation and annihilation operators are zero. 
Then, to the first order in P and to the g-th order in k 2 . the expectation value of \1Z\ 2 will be evaluated 

as 



1*1 



4iU / drjf I Im [u* (77' ')*]• P(r/)Re 



- %)Im KM*] • W)Re K(r/)*] 



(?) 



0(P 2 ), 



with 



*= [W(n,«o)] T 



(4.24) 



(4.25) 



Here we should remind that II is a pure imaginary number. 

In Sec. II V Al in order to obtain the solution solved in the forward direction using general slow-roll ap- 
proximation, we have assumed that all components of the multi-component scalar field are nearly massless. 
However, in Eq. I|4.24[) the first order term in P takes the form of Im [uq^)*] • P(7y)Rc [uq^)*]. Therefore, 
in this formula there is no need to assume all components are nearly massless. We only need to assume that 
the component of P connected to * is small. This seems a natural result, because basically the advantage of 
Wronskian formulation is that we only need to solve a single mode, as we mentioned before. We will prove 
the validity of this argument rigorously in our future paper [l^ . 



V. SINGLE FIELD 



Here, we consider the single scalar field case. In this case, we can easily obtain more explicit solutions 
described in an expansion with respect to fc 2 and g, which is correspond to P in Sec. IIVI respectively. 
Furthermore, it is known that the usual approach that uses 1Z C as a perturbation variable apparently breaks 
down when 4>' = 0. At that point special treatment is necessary to justify the evaluation of perturbation 
by means of 1Z C [13] • What we evaluate here is a similar quantity H c {r)f) = W but its argument is slightly 
different from lZ c (r]). Unless we choose rjf so that (j)'(r]f) = 0, lZ c (r]f) must be free from any singular behavior. 
This fact will become manifest below. 

A. Extended general slow-roll formula in the single field case 

In the case of single field, the equation for 5(f) f l|2.8|l reduces to 

{d 2 n + Q(r,) + k 2 )v(r ] ) = ^<p(r ] ), (5.1) 

with 

Q{ri) = ~, g( r] )=if--2, (5.2) 



9 



where z = a<fi' /Ti. 

The solution in the long wavelength approximation mS 1 ' satisfies 



o, a 



-k 2 n^ . 



The boundary conditions i|3.4[) at rj = 77 en d reduce to 

n = 0, n = — 1/z , 

where r] en d is taken as a time when scalar dominant phase ends. 
To the first order in k 2 , we have 

n« = n<°> + k 2 An^ 
r> drf 



(5.3) 



(5.4) 



-z{rf) 



- + fc^(r?) / rfT? 



loud 



Z 2 (t/) 



drf' 



As for the ^-expansion, now the lowest order mode function mq, which is a solution of 



is explicitly given by 



M 



//v - 



>7 



-ikr\ 



(5.5) 



(5.6) 



(5.7) 



Here we have chosen uq such that it corresponds to the positive frequency function of the Bunch-Davies 
vacuum, i.e. 



1 



u 



2k 



e~ ikr] for 



T] — > — oo. 



(5.8) 



In the present single field case, using Eqs. I|5.5|) and (|5.7(l . the power spectrum of 1Z C , which is given in 
Eq. I|4.24|) . to the first order in g and k 2 is more explicitly written as 



2k 2 af3) 7 



2 

X ~3 



drj — —T {rj ) 
if 



+2[a 1 + k 2 {aa + (i 1 )\ nk f" d^^Gtf) , 



(5.9) 



where 



a{rj) = - 



7fa) = —5- -a" 



1 




r] 2 




rf 


("4 


3 





m = - 



rj 2 I n K 



(o)\' i 



77A77, 



v 8 / n (0)\' ^4 / An (l)\' 



47(77) = — — r - — 
W7 30 V V 4 ' 



F{rj) = Re[F(77)] - % - 77 fe )Rc[F(7y)] (1) , g(n) = Re[G(r,)} - % - 7 7fc )Re[G(j 7 )] (1) 



(5.10) 



and where 



3i 2 
— (m Uo + Uqu ) = 1 + -/c 2 r/ 2 + 
?7 5 



G(t/) = &V (u Uo - u u ) = 1 + h*rf + 



(5.11) 
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Then, the asymptotic behavior of T(rj) and Q(rj) is 



lim F(rj) = - — sin(— krj) cos(— krf), 

kr/^-oo KT] 

lim Q(rj) = (— krj) 2 sin 2 (— krf), 

krj — > — oo 



lim J-(rj) 



35 



(-krj) 4 + ((-krj) 6 ), 



UmSfa) = -^(-kvf + 0((-krjf 



(5.12) 
(5.13) 

(5.14) 

(5.15) 



Here we show the plots of the window functions J-(rf) and G(ri)/rj 3 in FIG. ^ FIG. [I] tells that we have 
only to consider the contribution from the neighborhood of the horizon crossing time. Moreover, from these 
window functions' behavior, we can easily find that we may be allowed to set 77* to even if fall-off of g(rf) 
in the limit 77 — > is not very fast, as we mentioned in Sec. IIV Gl 



T( V ) 




- lo gin(~ k 'l) 



(kriy^r,) 











1 






J . , 




-2 


-1 
-1 


1 

- log 10 (- kri) 



FIG. 1: The window functions T(rj) and (krj) 3 G(rj) as functions of — log(— krj). 



B. 



When </>' = 0, it seems that the solution n' 1 ' diverges because z vanishes, but we will see below that in 
fact it does not. 

Using integration by parts, we can rewrite the integral / drj' / z 2 (rj') which appears in the Eq. 1)5.5(1 as 



J >7<md 



1 



-z(rj) 
1 



1 



z'(rj')z(rj' 
z(ij) 



z'(rj) z{rj en d)z'(rj end ) 



z'(rj') drj 1 \z(rj') 

n p 

-z(rj) I drj 1 

Vend J Vend 

+ z(rj) / drj' 



1 



1 d 
z(rj') drj' \z'(rj') 
z"(rj') 



z(n')z' 2 (rj') 



In a similar way, we also have 



An<°> = z(rj) f drj'^p- 
1 z(rj) 



z '(v) z(rj cnd )z'(rj cnd ) 

Jv.nd Z (V )Z 2 (V ) J n 



drj"z(rj")n m (v") + z(v) / 



(5.16) 



(5.17) 
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Here, z' and z" J z can be rewritten explicitly as 

/3 




11 6' 2 



4 H 2 m 4 



(5.18) 



Hence, we have z' ^ and z" / z ^ oo when <j>' = 0, as long as ^ 0. 2 Therefore we find that rjW doesn't 
diverge even if <j)' vanishes. 

Finally, it will be worth pointing out that the general slow-roll condition is even weaker than those shown 
in Ref. As noted earlier, under the general slow-roll condition, the slow-roll parameters, 6 = (f)/H<p : 
and 8/H can be the same order but both are assumed to be small in Ref. 0. However, we can rewrite 
g = 7] 2 z" I ' z — 2 by using the expression for z" j z given in Eq. (|5 . 1 8|) and 

rf = ^(l + 2e + 0(e 2 )), 

(5.19) 

where e = —H/H 2 . General slow roll assumes |<?| <C 1. When the background scalar field motion stops 
(4> ~ 0), the slow-roll parameter 5 diverges. Therefore, slow-roll conditions are violated at that point. 
However, from Eq. (|5.18|l . we can easily find that the condition \g\ <C 1 is maintained as far as V^/H 2 <C 1 
even when (f> ~ 0, namely, <j) ~ does not imply the breakdown of the general slow-roll condition. 

VI. SUMMARY 

We have proposed a new method for a systematic derivation of formulas for the spectrum of the curvature 
perturbations from multi-component inflation. First we have noted that the Wronskian W(n, f) between 
two solutions n and ip of the scalar field perturbation equation in flat slicing stays constant in time. Using 
this fact, we have shown that the evaluation of the curvature perturbations at the end of inflation reduces 
to an easier problem of solving a single "decaying" mode n backward in time from the end of inflation. 

We have also shown that, for any given perturbation <p the Wronskian W(n, <p) becomes identical to 
SN in the long wavelength limit by choosing an appropriate boundary condition for (n, n'), where SN is 
the perturbation in the e-folding number along the phase space trajectory of the homogeneous universe. 
More precisely, SN is a perturbation in the e-folding number from a given point in the phase space for the 
homogenous universe to a point when the trajectories converge to a unique one, with (ip,ip') identified as 
the perturbation in the background phase space of spatially homogeneous fields. 

Hence, with the aid of the SN formalism |l7lllSj . even if the phase space trajectories converge only after the 
scalar dominant phase ends, we can still evaluate the final value of the curvature perturbation by evaluating 
the Wronskian. Yet an open question is whether the perturbation in the true e-folding number, SAf, is well 
approximated by SN, which is determined by the dynamics of the homogeneous universe for the case with 
general matter fields. We plan to discuss this issue in our future publication 

As a good example to show the efficiency of the new formulation using the Wronskian, we have presented 
a quick derivation of an improved general slow-roll formula, whose original form was obtained in Ref. . In 
Ref. or in the general slow-roll formulae given in Ref. 1], it was assumed that super-horizon effects are 
negligible. However, it was pointed out in Ref. that there are cases in which they may not be negligible 
at all. 

The improvement that has been incoorporatcd for the first time in this paper is to use a systematic 
long-wavelength expansion for the late time evolution after horizon crossing. The matching of the two 
approximations, the general slow-roll approximation and the long wavelength expansion, has been easily 
performed by using the constancy of the Wronskian. We have shown that the formula which contains 



We can easily consider such cases, for example, the oscillating inflation proposed by Damour and Mukhanov [23 
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higher order corrections to the general slow-roll approximation can be obtained systematically. We have also 
presented an explicit power spectrum formula for the single field case. Furthermore, we have shown that 
our formula is valid even if the scalar field passes the point </> = where one of the slow-roll parameters. 6, 
diverges (6 = 4>/H<j) — > oo). This fact that any special treatment is not required even if <f> = occurs is a 
small but important advantage of our formulation in comparison over previous treatments ■ 
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